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ABSTRACT 

We propose that the curious Galactic Center filament known as "The Snake" 
is a twisted giant magnetic flux tube, anchored in rotating molecular clouds. 
The MHD kink instability generates coils in the tube and subsequent magnetic 
reconnection injects relativistic electrons. Electrons diffuse away from a coil at an 
energy-dependent rate producing a flat spectral index at large distances from it. 
Our fit to the data of Gray et al. (1995) shows that the magnetic field ~ 0.4mG is 
large compared to the ambient ~ 7/iG field, indicating that the flux tube is force- 
free. If the relative level of turbulence in the Snake and the general interstellar 
medium are similar, then electrons have been diffusing in the Snake for about 
3 x 10 5 yr, comparable to the timescale at which magnetic energy is annihilated 
in the major kink. Estimates of the magnetic field in the G359. 19-0.05 molecular 
complex are similar to our estimate of the magnetic field in the Snake suggesting a 
strong connection between the physics of the anchoring molecular regions and the 
Snake. We suggest that the physical processes considered here may be relevant 
to many of the radio filaments near the Galactic Center. We also suggest further 
observations of the Snake and other filaments that would be useful for obtaining 
further insights into the physics of these objects. 

Subject headings: Galaxy: center, ISM: kinematics and dynamics, ISM: magnetic 
fields, magnetohydrodynamics, radio continuum: ISM, stars: formation 

1. Introduction 

The fundamental nature of the numerous radio filaments observed near the Galactic 
Center (see LaRosa et al. (2000) for an overview) is unclear. In particular: How are they 
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formed? What are the sources of the relativistic electrons? Why do they have such large 
magnetic field strengths? What is the reason for the unusually flat spectral indices in 
some of them? The discovery of the Snake (Gray et al. 1991), a 60 pc by 0.4 pc filament, 
approximately 150 pc to the West of Sgr A, resulted in the additional feature that the 
characteristics of the radio emission are strongly related to a morphological "kink" . Further 
observations (Gray et al. 1995), revealed a major and a minor kink, with the radio intensity 
and spectral index systematically varying away from the major kink. Various models have 
been proposed to explain the spectacular elongated structure of the Snake (see Gray et al. 
(1995)). However, none of them have been able to successfully explain its major features. 
In our view, the Snake is in many respects one of the least complicated of the filamentary 
features close to the Galactic Center and a good physical model may hold the key to 
understanding the numerous arcs and filaments in the Galactic Center region. 

We propose that the Snake is a magnetic flux tube with both ends anchored in dense 
rotating material (molecular clouds and/or associated HII regions). Differential rotation of 
the tube at either or both ends produces a monotonically increasing toroidal magnetic field 
and when this reaches a critical value "coils" or "loops" are formed through localized kink 
instabilities. Release of the magnetic energy stored in a coil, through magnetic reconnection 
(and possibly shocks), is the source of the relativistic electron energy. Energetic electrons 
diffuse away from each coil at an energy- dependent rate causing a flattening spectral index. 

2. Details of the model 

Let us take R to be the radius of the tube, L tu be its length and and B z the toroidal 
and axial components of the magnetic flux density, with magnitude B, in a cylindrical 
(r, 0, z) coordinate system coincident with the central axis of the unperturbed tube. As the 
end(s) of the tube is(are) rotated, the toroidal field increases and when B^jB z ~ 1, the 
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kink instability produces a coil of magnetic flux with magnetic energy AE rn A^ttR^B 2 
(Alfven (1950), pll7). We identify the observed major kink with such a coil; the minor 
kink may be another coil at a different stage of development. Numerical simulations, in a 
solar physics context (Bazdenkov & Sato 1998; Amo et al. 1995), have shown that the coil's 
magnetic energy is annihilated on a time scale of order a few transit Alfven crossing times, 
L tuhc /vA, where va is the Alfven speed. Magnetic reconnection and possibly associated 
shocks in the coil can provide a source of energy for the acceleration of electrons to 
relativistic radio-emitting energies 2 . Once accelerated, the electrons diffuse and we model 
the diffusion of particles away from the acceleration region by a one-dimensional diffusion 
equation, assuming a uniform cross-section for the tube. Taking f{p) to be the electron 
phase-space density, x the spatial distance along the flux tube away from the major kink, 
K{p) the momentum-dependent diffusion coefficient, and C(p) the rate per unit volume of 
momentum space at which electrons are injected into the coil by the acceleration process, 
we adopt the diffusion equation: 

df(p,x,t) 9r {p) df( P x,t) \ cmx) (1) 



dt dx \ dx 

The delta function indicates that we are treating the coil, located at x — 0, as an 
infinitesimally small volume with respect to the rest of the tube. We also assume that 
the timescale for dissipation of energy in the coil is greater than the time over which the 
electrons diffuse. We discuss this assumption further below. We adopt a power-law for 
K(p), i.e. K{p) = Kolp/po) 13 , where K is the diffusion coefficient at an electron momentum, 
p of a GeV/c. We also assume that C{p) = C (p/p )~ s . Integrating equation (1) from 
x = 0~ to x = + gives the boundary condition at x — 0, namely, 

df(p,x,t) lC(p) 



dx 



x=o 2 K{p) 



(2) 



2 The presence of shocks cannot be inferred from the Bazdenkov et al. simulations since 
they assumed an incompressible fluid. 
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We normalize the diffusion equation and the boundary condition by introducing the 
distance L = 1 5 pc between the major and minor kinks as a fiducial length together 
with the normalized variables r = K(p)t/L 2 = (K t/L 2 ) (p/po)^, £ = x/L and 
g(£,,T) = 2K C 1 L~ 1 (p/po) s+l3 f(p,x,t). In these variables, the diffusion equation and the 
boundary condition become 



dg d 2 g 



with 



dg_ 



5=0 



(3) 



for which the solution is 



1/2 



1 - erf 



(4) 



The number density of particles per unit Lorentz factor, 7, is 



N( 1 ,Z,T)=N ol - a g(Z,T), 



(5) 



where iV = 27r(m e c) 3 (CoL/fT )7o +< \ 7o = (m e c)"Vo and a = s + (3 - 2. We use ^(7) to 
evaluate the angle- averaged synchrotron emissivity, < j u >, as a function of frequency, u, 
from 



< >= (V3e 2 /4c) N v B (2v/3v B f- a)/2 / y^ 2 g(^ r)F{y) dy 
Here, v B = eB/(27rm e c) and the angle-averaged single electron emissivity 
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F(y) — V a/I — y 2 /u 2 if 5/3 (u) du. Because of the 7 dependence in <?(£,r), the spectral 
index near the kink is not simply (a - l)/2 but (s + (5/2 - 3)/2 = (a - l)/2 - /3/4. 

For a uniform flux tube, of radius R and angular diameter imaged with a circular 
Gaussian beam with standard deviation, a, the angle-averaged flux density per beam is 
F v = 4(2tt/3) 1 / 2 ((7$) [I ($ 2 /12a 2 ) + Ji($ 2 /f2 ( j 2 )] (< j v > R) where J and h are modified 
Bessel functions. Using this expression and the expression for the emissivity, we have 
performed a least squares fit to the observed flux densities (Gray et al. 1995) at 1446 and 
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4790 MHz, for which, a = 4.90", R = 0.22 pc and $ = 9.4" are appropriate. We restrict 
the fit to the region of the Snake between galactic latitudes of —3' and —17' that is clearly 
related to the major kink; the radio images and the spectral index plot (Gray et al. (1995), 
Figure 12) indicate that the minor kink influences the flux density south of b = —17'. The 
parameters of the fit (see Figure 1) are N = 3.5 x 10~ 5 cm -3 , B = 0.37 mG, a = 2.14, 
(3 = 0.57 and r = K t/L 2 = 0.46. The deviations of the fit from the model are most likely 
the result of systematic effects, such as variation in the field strength and width of the tube, 
rather than the statistical uncertainty in the data. The magnetic field can be estimated 
from the data, since it is related to the frequency index, which varies significantly along 
the tube. However, experimentation with the fit shows that the precision of the estimate 
of B is not much better than a factor of a few although the result that the magnetic field 
is much greater than the ambient interstellar value ps 7/iG (Gray et al. 1995), is robust. 
Therefore, the flux tube is force-free. 

If we assume that the injected electron spectrum extends between momenta, p\ and 
P2, the total electron energy in the flux tube can be estimated as follows. Let A be the 
tube's cross-sectional area, then the rate, P e of relativistic electron energy injection into the 
tube, is P e = AttAc J^p 3 C(p) dp = 4irAC cp 4 (4 - s)" 1 [(p 2 /po)^ s - (piM)) 4 ^] The total 
relativistic electron energy in the flux tube is therefore, 

E c , tot = P c t = 27 ri? 2 m e c 2 7 2 ^ (fr) (^) [^/p^ s - [p./p^] 

« 6.3 x 10 44 (p 2 /po)° A3 ergs (7) 

using the expression for iV (following equation (5)) and the derived parameters. Time 
enters through the dimensionless model parameter, r = K t/L 2 . For an upper cutoff of 
10 GeV, £ e ,tot ~ 3.9 x 10 45 ergs and for 100 GeV, E tot ^ 1.1 x 10 46 ergs. The available 
magnetic energy stored in the coil, Eb ~ 3.5 x 10 46 ergs, exceeds the total energy in 
electrons, as required for consistency of the model, but not necessarily by a large factor, 
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depending upon the maximum energy to which electrons are accelerated. This is a 
consequence of the high value of the derived magnetic field. These comparisons also suggest 
that the diffusion of electrons has been taking place for a time comparable to the total time 
available to annihilate the coil, consistent with our not observing the process at a special 
epoch and consistent with the assumption of continuous injection. 

The energy density of electrons, e c , can be derived simply from equation (5) for Ni^y). 
Near the coil (f = 0), 

ee = _^]t + 2 i^/P^ a+m+2 - ^/Por a+m+2 ] (8) 

where N x = 2tt' 1 / 2 N Tq /2 1o P/2 . The energy density is insensitive to the upper cutoff, and 
for p 2 = 10 GeV/c, e c ~ 6 x 10~ n ergs cm~ 3 , an order of magnitude larger than the energy 
density of the ISM. 



3. The spectrum of hydromagnetic turbulence 

It is usually assumed that relativistic electrons resonantly scatter off a pre- 
existing level of hydromagnetic turbulence since resonant Alfven waves are damped 
rapidly in the warm interstellar medium by ion-neutral collisions, although this 
constraint is more important at higher then GeV energies (e.g. Melrose (1982)). 
Let the energy density per unit wave number, k, of resonant Alfven waves (either 
pre-existing or self-excited) be W(k) = Woik/ko)^, where k = eB/cpo corresponds 
to waves resonating with 1 GeV electrons. Then, the spatial diffusion coefficient 
for relativistic electrons is K(p) = 497(77 + 2)/9tc (c 2 /e) B~ 1 p[W m /k R W(k R )] where 
W m = B 2 /8n and k R = eB/cp is the resonant wave number (Melrose 1982). Numerically, 
K(p) = 1.4 x 10 19 r7(r7 + 2) (5/mG)" 1 (W m /k W ) (p/pof" cm 2 s" 1 . Hence, (3 = 2 - r] and 
i] ?s 1.43 for the parameters of our model. Our value of j3 ~ 0.57 is close to the Ormes & 
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Protheroe (1983) value of 0.7 derived from a cosmic ray propagation model. However, more 
recent models adopt a Kolmogorov value, j3 = 1/3, combined with the effects of "minimal 
reacceleration" (e.g. Ptuskin et al. (1999)). We have yet to take account of this possible 
effect. In order to constrain f3 more effectively, one needs to take into account, not only 
minimal reacceleration, but the other physical parameters in the problem such as the width 
and strength of the flux tube and time variability in the injection process. 

The time since the coil started to inject electrons along the flux tube is of interest for 
comparison with other timescales. This is t — L 2 /K r « 3.1 x 10 5 K~\ % (ro/0.46) yr. For 
cosmic rays, Ptuskin et al. (1999) take K ~ 10 28 cm 2 s -1 , implying that k Q W Q /W m ~ 10~ 6 . 
If the same relative level of turbulence exists in the Snake, then K ~ 10 26 cm 2 s" 1 since 
K oc B^ 1 and t ~ 3 x 10 5 yr. According to the numerical simulations of comparatively 
short wound flux tubes (Bazdenkov & Sato 1998), a coil disappears "explosively" in 1 — 3 
Alfven-crossing times, L tuhc /vA- For the Snake, this is 1.5 — 4.5 x 10 5 (n/10 cur 3 ) 1 / 2 yr 
based on L tubc pa 60 pc and a number density n ~ 10 cm -3 for the ambient interstellar 
medium (Gray et al. 1995). Explosive bursts recur approximately every 5 Alfven times, i.e. 
approximately every 7.5 x 10 5 (n/10 cm -3 ) 1 / 2 yr. Thus, for K ~ 10 26 cm 2 s -1 , the time over 
which the electrons have been diffusing is of order the timescales of the energy releasing 
process and, as above, we are not required to invoke a special epoch of observation. This 
also again justifies our assumption of continuous injection over the diffusion timescale. 



4. The origin of the magnetic field 

We suggest that a twisted magnetic flux tube with the properties required by our 
model would arise in the following way: The magnetic flux tube is initially anchored at both 
ends in molecular clouds before they undergo contraction and initiate star formation. Since 
molecular clouds condense from the warm interstellar medium, the cloud magnetic field at 
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the pre-contraction phase exceeds the general ISM value. Hence such a flux tube emerging 
into the ISM from a molecular cloud would expand. As one or both of the clouds contract 
and rotate, a significant toroidal field is produced and the resulting magnetic curvature 
force draws in the flux tube, thereby increasing the poloidal flux density in the entire flux 
tube to the molecular cloud value and eventually leading to a force free state which is also 
unstable. (Both a force-free configuration and instability require at the boundary of the 
flux tube to be of order B z .) 



There is reasonably good evidence for the anchoring of the Snake in molecular clouds 
or associated HII regions. Observations by Uchida et al. (1996) reveal that the northern 
end of the Snake intersects an HII region in a CO "hole" in the molecular cloud and HII 
region complex, G359. 19-0.05. The supernova remnant observed in projection against the 
southern end (Gray et al. 1995) provides at least circumstantial evidence for molecular 
material in that region. Another possibility is that the Snake could be part of a giant loop, 
the other end of which is anchored in another region of the Galactic Center. The generation 
of toroidal field is inextricably linked to the evolution of magnetic field and angular velocity 
in contracting molecular clouds. Here we discuss some of the physics involved. However, 
it will be evident that a complete analysis involves substantial issues in the physics of 
contracting molecular clouds that are beyond the scope of this letter. 

The generation of toroidal fields by contracting clouds is related to the radiation of 
torsional Alfven waves by the rotating field anchored in the cloud (see Mestel (1999), p 
453). Let p be the mass density of the background ISM and let B z be the (uniform) 
poloidal flux density of the cloud and ISM. The angular velocity, f2 and toroidal component, 
B^ of the field in the tube are governed by: 

on b z ob^ ob^ D on 
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The angular velocity in the background medium in which the Alfven speed is va, satisfies 

W = ^a* (10) 

As a consequence of these equations, the toroidal field in a tube which is anchored only in 
the cloud is « -r(47rp ) 1/2 fio ~ -1-5 x 10- 6 (n/10 cm- 3 ) 1 /2(Q /km s" 1 pc" 1 ) G, where 
flo is the angular velocity of the cloud (Mestel 1999). In order to generate ~ 4 x 10~ 4 G, 
one requires an extraordinarily large value of fl ~ 270 km s _1 pc _1 . 

Consider now a tube anchored, as well, in another molecular cloud. The boundary 
conditions on Q are now, Q = Q at z = and Q = at z = L tu bc- On a timescale long 
compared to the Alfven time, the relevant solution of equation (10) is Q = Q (l — Z /L tu bc), 
implying from the second of equations (9), that 

dB <i> _ r> m 

— --B - (11) 

Ol i^tube 

Further, if Q is constant for a time t, then, after the flux tube has been twisted through 
an angle A0 = f2 t we have = —rA(pB z /L tu ^ c . This expression can be derived 
simply from flux freezing in a twisted flux tube, neglecting the propagation time of 
Alfven waves and is used by Alfven (1950), for example, in the theory of the instability 
of twisted flux tubes referred to above. The purpose of the derivation here is to make a 
clear connection with the torsional Alfven waves involved in the standard treatment of 
magnetic braking. On the basis of this model, ~ B z would be attained for the Snake if 
ilt ~ L tuhc /R PS 300(i?/0.2 pc) _1 radians. This would be the case, for example, for a cloud 
rotating at an angular velocity of 30 km s _1 pc _1 for a period of 10 7 yr. 

However, the above estimates of the toroidal field are somewhat contrived in that they 
assume that the flux tube has a constant radius which is clearly not the case in a contracting 
cloud. Moreover, the poloidal field changes under the opposing effects of compression and 
ambipolar diffusion and the evolution of the field in the external ISM (see above) is also 
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important. However, the main point is that the anchoring of the magnetic field at both ends 
substantially affects the generation of toroidal field and the requirements on the angular 
velocity. 

A related issue involves that of subcritical or supercritical contraction of the anchoring 
cloud(s). In subcritical contraction, the initial magnetic field supports the cloud and 
contraction occurs as a result of ambipolar diffusion of the magnetic flux together with 
magnetic braking. The latter prevents the cloud from approaching centrifugal equilibrium 
(e.g. Basu & Mouschovias (1995) and references therein) and the angular velocity remains 
low, except for the latest, rapid stages of contraction when the core becomes critical. On 
the other hand, in the supercritical case, the magnetic field cannot halt contraction, and 
if the cloud has an initial angular velocity, angular momentum conservation causes it to 
spin up, until it reaches centrifugal equilibrium. Further contraction is then the result of 
magnetic braking (Mestel & Paris 1984). Therefore, supercritical contraction may offer 
the best prospect for twisting of the flux tube. However, this is a complex issue and 
investigation is deferred for future work. Note also, that Basu & Mouschovias (1995) began 
their simulations with slowly rotating clouds. 

An independent estimate of the poloidal magnetic field is of interest since this is one 
of the key parameters in our model and also determines whether the molecular cloud 
contraction associated with the Snake is sub- or supercritical. We use parameters derived 
by Uchida et al. (1996) who estimate that the clouds adjacent to the CO hole have masses 
~ 5 x 1O 3 M with radii, R ~ 2 pc. If we assume that these clouds have formed through 
compression of spheroidal regions (with semi-axes _R and KR ) of the warm interstellar 
medium of density n ~ 10 cm" 3 and with an ambient magnetic field of B ~ 7 x 1CT 6 G 
(Gray et al. 1995), then flux conservation implies that 
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B z 
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R? \_4nK /xn m p 
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" 7 * - 4 (rT^c) (^P^" (s^f <™ 

Using expressions derived in Mestel (1999) (pp. 429 et seq.), the critical magnetic field 

B ciit « 6 x 1(T 5 - 1 x 1(T 4 f ^ ( G (13) 

\5 x 1O 3 M ; V 2 Pc/ V ; 

with the lowest (highest) value for a spherical (flattened) cloud. These estimates of B z and 

B CTit raise the prospect that the clouds in the G359.19-0.05 complex are subcritical unless 

the original ISM region is highly prolate (K » 1). However, our main point is that these 

indicative estimates of both the actual and critical fields are similar to our estimate of 

4 x 1CT 4 G from our model for the Snake. It is therefore possible that the flux density in 

the associated contracting molecular cloud is not greatly enhanced over its initial value, as 

a result of ambipolar diffusion. Many of the molecular clouds in the vicinity may have fields 

of this magnitude, but it is not until they are twisted and "lit up" by reconnection induced 

by an instability that they become observable. Interestingly the Basu & Mouschovias 

(1995) simulations show that the magnetic field outside the supercritical core is indeed time 

independent. 

5. Discussion 

High resolution images of the Snake near the major kink (Gray et al. (1995), Figures 
10 and 11) show that the Snake appears to be split in two, similar to flux tubes in the 
late stages of the Bazdenkov & Sato (1998) simulations. This provides additional support 
for our proposal. As well as being directly applicable to the Snake, our model may open 
up a number of appealing possibilities for the dynamics of other magnetic filaments in the 
interstellar medium near the Galactic Center. It is feasible, in an environment with such a 
large density of molecular clouds, that rotation of magnetic flux tubes threading their cores 
would lead to kink instabilities and reconnection. The flattening spectral index of the Snake 
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depends upon the differential diffusion of energetic particles and this is such a well known 
phenomenon in the regular ISM that its role in this model is unremarkable. The index of 
the momentum dependence of the diffusion parameter that we have inferred for the Snake 
is in the range of values used in models of cosmic ray propagation. More detailed models 
raise the prospect of a better determination of this index as well as the other parameters in 
this model. In particular, a more precise value of the index of the momentum dependence 
of the particle creation rate, which our current estimates place near the value associated 
with strong shocks, would be of interest. The type of data that is essential for more detailed 
modeling include (1) High resolution images showing more clearly the structure of the flux 
tube over its entire length. (2) Flux densities at different frequencies, enabling one to better 
constrain the energy dependence of the diffusion parameter and the magnetic field. 

It remains to be shown that flux tubes in the Galactic Center can be wound up to such 
an extent that they produce a significant toroidal field. If an anchoring cloud acquires an 
angular velocity which persists for long enough that it rotates through 300 radians (~ 50 
revolutions), that would be sufficient. Another possibility is that the required winding 
could be produced in a contracting cloud, although this idea requires further investigation. 
In our independent estimates of the magnetic field in the Snake and in the related issue of 
subcritical or supercritical contraction we have concentrated exclusively on the molecular 
complex at the northern end of the Snake. Currently, however, there is little known about 
the southern end and the masses and rotation rates of molecular clouds in that region will 
certainly be of interest. 

We thank Dr. R. Protheroe for helpful information on cosmic ray propagation, Prof. L. 
Mestel for advice concerning the physics of molecular clouds and the referee for constructive 
comments on the original version of this paper. 
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Figure Caption 

Figure 1: The model fit to the 1446 MHz (filled squares) and 4790 MHz (open circles) 
flux densities (Gray et al. 1995) for the Snake. The model fit only applies to the region 
associated with the major kink. 
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